Driving Bose-Einstein condensate vorticity with a rotating normal cloud 
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We have developed an evaporative cooling technique that accelerates the circulation of an ultra-cold 
8r Rb gas, confined in a static harmonic potential. As a normal gas is evaporatively spun up and 
cooled below quantum degeneracy, it is found to nucleate vorticity in a Bose-Einstein condensate. 
Measurements of the condensate's aspect ratio and surface-wave excitations are consistent with 
effective rigid-body rotation. Rotation rates of up to 94% of the centrifugal limit are inferred. A 
threshold in the normal cloud's rotation is observed for the intrinsic nucleation of the first vortex. 
The threshold value lies below the prediction for a nucleation mechanism involving the excitation 
of surface-waves of the ground-state condensate. 
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To paraphrase an ancient riddle, what happens when 
an irresistible torque meets an irrotational fluid? The 
answer has been known for more than 50 years: a quan- 
tized vortex is nucleated. Vortices alone contribute to a 
superfluid's rotation, so that the bulk of the fluid may 
remain curl-free. The nucleation of vortices in bulk su- 
perfluid Helium has been the topic of extensive study 
(for a review see m). In the archetypical experiment, 
a rotatable pot filled with a mixture of superfluid and 
normal liquid Helium undergoes gradual angular acceler- 
ation. The normal fluid and the walls of the pot rotate 
together as a rigid body, defining a rotating environment. 
At some threshold angular velocity, a vortex line is nu- 
cleated at the circumference of the pot, and then quickly 
migrates inward until it is collinear with the axis of ro- 
tation. Further angular acceleration results in the nu- 
cleation of more vortices; eventually the fluid is filled 
with an array of vortex lines A central theme jj| 
of this research is the question: to what extent is the 
nucleation process "extrinsic," i.e. dependent on such 
details as the roughness of the surface of the walls, and 
to what extent is it "intrinsic" i.e. driven (in the 

limit of microscopically smooth walls) by the flow of nor- 
mal fluid along the boundary of the superfluid? In the 
analogous rotating-potential experiments with a dilute- 
gas Bose-Einstein condensate (BEC), the confining po- 
tential and the normal fluid typically rotate at different 
rates In this context, the extrinsic-intrinsic question 
can be restated as: is it the confining potential or the 
normal fluid that defines the rotating environment? 

Vortices in a BEC have been created with wavefunction 
engineering Q, through the decay of solitons and 
in the wake of moving objects p^fO] . The first rotating- 
potential experiment to detect vortices in a BEC was 
performed by the Paris group |l2| ; results have also been 
obtained by the MIT Q and Oxford @ groups. In 
these experiments the role of the normal fluid was sec- 
ondary to that of the rotating potential; it is conceiv- 
able the normal fluid was not rotating at all. This paper 



presents vortex nucleation experiments performed in the 
opposite limit, namely in the environment of a rotating 
normal gas in a static confining potential. Such an en- 
vironment allows for the isolated study of the intrinsic 
mechanism for vortex nucleation. 

Our experiments begin with a magnetically trapped 
cloud of about 6 • 10 6 87 Rb atoms, in the \F — 1, mp = 
— 1) hyperfine state, cooled close to the critical temper- 
ature T c = 67 nK. The atoms are initially confined in 
an axially symmetric, oblate and harmonic potential |l5| ] 
with axis of symmetry along the vertical ('z') axis. To 
induce rotation of the cloud, we first gradually apply 
an elliptical deformation to the potential in its horizon- 
tal plane of symmetry, and then rotate the deformation 
Jl6| about the vertical axis at a fixed angular frequency. 
The rotating potential is characterized by an axial fre- 
quency lu z — 27r(13.6)Hz, time-averaged radial frequency 
(u) p ) = 27r(6.8)Hz and a horizontal ellipticity of 25%. 
Such a large rotating trap asymmetry, accessible in the 
oblate configuration of our apparatus, is found to be 
necessary not only to get the cloud rotating, but also 
to sustain ongoing rotation. Moreover, in steady state 
the cloud does not reach the rotation rate of the applied 
asymmetry. We believe that the stirring process is fight- 
ing a small, static asymmetry that acts to despin the 
cloud 0]. 

In thermal equilibrium, a normal cloud rotates as a 
rigid body with the centrifugal force causing the cloud 
to bulge outwards in the radial direction. In order to de- 
tect the rotation, we use a nondestructive phase-contrast 
technique to image the cloud in situ from the side. Four 
sequential pictures of a given cloud are taken to average 
over oscillations in the widths and to improve the signal 
to noise. The cloud temperature is extracted from its ver- 
tical width a z , which is unaffected by rotation about the 
vertical axis. The rotation of the cloud £In is determined 
from the aspect ratio, A = o~ z /a Pl using the relation: 

n/ujp = v 7 ! - (A/A ) 2 (1) 
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where A D is the static aspect ratio. The technique of 
side-view imaging is crucial for distinguishing between 
changes in radial size due to temperature and to rotation. 

For stirring rates up to 2.5 Hz, the rotation of the 
cloud reaches its steady-state value by 15 s or less. After 
15 s, the rotating trap asymmetry is ramped off, leaving 
non-condensed clouds a factor of 1.2 - 1.3 above T c , for 
stirring frequencies - 2.5 Hz. Radio frequency (rf) evap- 
oration is then used to cool the normal cloud to BEC. In 
the oblate configuration of our TOP trap |l5"| , we have 
found that rf evaporation immediately quenches the ro- 
tation, presumably because the selection process, which 
removes atoms with large radial displacements, preferen- 
tially removes atoms with large axial components of their 
angular momentum. By adiabatically distorting the trap 
into a prolate geometry with {u p , lj z } = 27r{8.35, 5.45}Hz 
JUIl , we can instead cool the cloud by removing atoms 
with large axial displacements and thereby reduce the 
effect of the evaporation on the axial angular momen- 
tum. As the normal cloud is evaporated, its aspect ratio 
is observed to decrease continuously, indicating a mono- 
tonically increasing rotation rate (Fig. |l]a). During the 
evaporation, the angular momentum per particle of the 
normal cloud remains roughly constant, even though the 
number of atoms is reduced by over a factor of five and 
temperature, by a factor of four (Fig. [j]b). As the cloud 
cools and shrinks, it must spin up for the angular mo- 
mentum per particle to remain fixed. 

To reach significant rotation rates by the end of evapo- 
ration requires the lifetime of the normal cloud's angular 
momentum to be comparable to the evaporation time. 
The nearly one-dimensional nature of the evaporation to- 
gether with the low average trap frequencies make cooling 
to BEC in the prolate trap very slow (~50 s) . We obtain 
angular momentum lifetimes this long by shimming the 
azimuthal trap symmetry to better than 0.1%. 

Towards the end of the evaporation, a condensate be- 
gins to appear at the center of the rotating normal cloud 
||19| . We discuss first the results of experiments in which 
we continue the evaporation until little or no normal frac- 
tion remains. In this case, we find that the rotating nor- 
mal component has given birth to a condensate distended 
in its radial dimension, as one would expect for a classical 
rotating body under the influence of the centrifugal force. 
This effect is reminiscent of liquid Helium experiments, 
in which the surface of a rotating bucket of superfluid 
exhibits the same meniscus curvature as for an ordinary 
viscous fluid p0| . For large enough numbers of vortices in 
the condensate, the correspondence principle would sug- 
gest that the rotation field, coarse grained over the cloud, 
should go over to the classical limit of rigid-body rota- 
tion. In this limit, the classical Eq. ([!]) should connect 
condensate aspect ratio to rotation rate. 

Alternatively, we can study the angular momentum in 
the BEC directly by exciting quadrupolar surface waves 
pT|,p2| with a rotating weak deformity of the magnetic 
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FIG. 1. (a)Aspect ratio of a rotating normal cloud during 
evaporation preferentially along the axis of rotation. Three 
different initial rotations of the cloud are shown, obtained by 
first stirring the cloud for 15s with an applied rotation rate 
of Oujp (inverted triangles), 0.07lo p (open circles) and 0.37oj p 
(squares), (b) Angular momentum per particle (squares) and 
temperature (triangles) of the normal cloud during evapora- 
tion for the 0.37u p case in (a). 

trap. The quadrupolar surfaces waves are characterized 
by angular momentum quantum number m z = ±2 where 
the m z = +(— )2 excitation is taken to be co- (counter- 
propagating with the rotation of the condensate. By 
varying the initial stir rate applied to the normal cloud, 
condensates of different aspect ratio can be accessed. In 
Fig. ||a, the frequency of the m z = ±2 modes is shown 
as a function of condensate aspect ratio. The m z = +2 
mode is seen to speed up and the m z — — 2 to slow down 
due to the presence of vorticity in the condensate. 

For small rotation rates, the splitting between the 
m z = ±2 modes is predicted to be linearly propor- 
tional to the mean angular momentum of the condensate 
p3| , p4j . In the large- SI limit of rigid-body rotation, Zam- 

have used a sum-rule argument 



belli and Stringari |24p; 
to show that the splitting between the modes is simply 
2f2, and, further, that the sum of the squared frequen- 
cies of the two modes is independent of rotation rate. A 
best fit of this model to the combined m = ±2 data is 
shown in Fig. ^a, where the rigid-body rotation rate has 
been inferred from the condensate aspect ratio and the 
classical Eq. (Q). Perhaps more intuitively, the frequency 
splitting is plotted explicitly versus inferred rotation rate 
in Fig. [|b. The excellent agreement with the model of 
Zambelli and Stringari is compelling evidence in favor of 
the reasonableness of using Eq. ([!]) to connect the con- 
densate aspect ratio with its effective rotation rate. This 
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FIG. 2. Quadmpolar surface-wave spectroscopy of con- 
densates formed in a rotating normal cloud. (a)Quadrupolar 
frequency as a function of condensate aspect ratio for the 
m = +2 (squares) and m — —2 (triangles) surface waves. 
Solid lines are a single fit to the combined data using the the- 
ory of Zambelli and Stringari J2^,^]. Dotted lines indicate 
average frequency and aspect ratio for a static BEC. (b)The 
splitting between m = ±2 frequencies scaled by twice the 
trap frequency, plotted explicitly as a function of BEC ro- 
tation rate inferred from the aspect ratio. The solid line is 
the prediction from the same theory as in (a). Each plotted 
point is obtained from a single m — —2 measurement in (a) 
combined with a spline interpolation to the relatively quiet 
m = +2 data. Only aspect ratios smaller than 1.43 (corre- 
sponding to Qbec I Up > 0.35) are included to avoid obtaining 
imaginary rotation frequencies due to experimental noise in 
the aspect ratio. 

is further born out by extensive 3-D numerical simula- 
tions of the Gross-Pitaevskii equation for the parameters 
of our experiment, by Feder and Clark. Their numerical 
simulations confirm that a condensate in an environment 
rotating at frequency fl > 0.5 lo p will equilibrate close to 
the aspect ratio given by Eq. (1) p6[ . 

In pure condensate samples we have observed aspect 
ratios as pronounced as 0.35 A Q , corresponding to a rota- 
tion rate of 0.94 u> p . The rapid rotation rate, combined 
with the increased condensate area arising from its radial 
bulge, mean that the condensate must be supporting a 
large number of vortices. Feder and Clark [^6| calculate 
56. With these initial conditions, we have observed con- 
tinued rotation for at least 140 s. 

If the evaporation is stopped before the normal cloud 
has been completely removed, a comparison can be made 
between the aspect ratios, and hence rotation rates, of 
the condensate and normal cloud. By adjusting the ini- 
tial stir rate applied to the normal cloud and the depth 
of the evaporation, we are able to reach different rotation 
rates for a given condensate fraction. After the evapo- 
ration is stopped, a time of 5 s is allowed for the gas to 
rethermalize to the last evaporative cut. We then take 
four nondestructive pictures and fit the images to a two- 
component distribution. Figure |^ shows a plot of the 
condensate aspect ratio Xbec compared with the aspect 
ratio A at of the normal cloud where each point represents 
a single realization of the experiment. 
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FIG. 3. (a) Aspect ratio of the BEC vs. that of the normal 
cloud after evaporation halted. Right and top scales provide 
a conversion from aspect ratio to classical rigid-body rotation 
rate for BEC and normal cloud respectively. Data for evap- 
oration of both a static (empty squares) and rotating cloud 
(filled squares) are shown. Dotted lines indicate average static 
aspect ratios for both BEC and normal cloud. A solid 1:1 line 
is superimposed on the data. Three representative integrated 
density profiles of two-component clouds indicate the range 
of different aspect ratios observed. (b)A magnified version 
of the region of high aspect ratio (low rotation rate) in (a). 
Added to the plot are data (triangles) obtained with evapo- 
rative spin-up 3 times slower than for the filled squares. A 
dashed line indicates the aspect ratio expected for a BEC with 
a single, centered vortex as calculated by Feder and Clark [M . 



3 



For \n < 1.36 (fijv > 0.5 uj p ), the condensate aspect 
ratio closely tracks the normal aspect ratio, providing 
a further manifestation of the correspondence principle 
for a highly rotating BEC. In the vicinity of Xn = 1-48 
= 0.35 Lo p ), we observe threshold behavior in the 
condensate rotation. At this low value of rotation we 
don't expect the rigid-body model to be valid for the 
condensate; but we make use of a numerical calculation 
by Feder and Clark |^6| to indicate the change in conden- 
sate aspect ratio associated with the presence of a single, 
centered vortex (Fig. ||b). There is considerable scatter 
in the data so one cannot make a strong statement about 
the nature of the threshold shape, but clearly, somewhere 
between 0.32 < Qn/uj p < 0.38 the first vortex is nucle- 
ated. 

A comparison of this threshold location with two the- 
oretical rotation rates for our particular cloud provides 
some insight into the nature of the vortex nucleation. 
The first value is the threshold for the thermodynamic 
stability of a single vortex. For our experiment this 
number is 0.2 - 0.25 lo p |2?|| , distinctly lower than our 
observed threshold for nucleation. The second value is 
Umim the frequency at which the slowest surface- wave 
mode propagates around the circumference of the con- 
densate. The Paris group has shown that for their "ex- 
trinsic" nucleation process (vortices nucleated by a rotat- 
ing asymmetric potential) the key mechanism is the non- 
linear excitation of surface waves jl2|,^8|. Results from 
MIT Jll],[l3| and fr om Oxford ]l4] | arc also consistent with 
such a mechanism. For the parameters of our experiment, 
u m in = 0.4 Lo p [ p9[ . Our observed threshold for "intrin- 
sic" nucleation is clearly under this value, and thus we 
cannot interpret our effect in terms of a normal "wind" 
exciting surface waves on our condensate. 

The data presented in Fig. [| include a range of con- 
densate fraction from 0.1-0.6 for each rotation rate of 
the normal cloud, although no segregation relative to ei- 
ther axis is evident for plots of different BEC fractions. 
Moreover, by reducing the rate of evaporation, we have 
decreased the rate of acceleration of the normal cloud ro- 
tation by a factor of 3, and still observe a threshold for 
vortex formation between 0.32 and 0.38 (Fig. ||b). 

In future work we plan to use the controlled rotation 
of the normal cloud to study the life cycle of a BEC's 
vorticity [j30| , including any hysteresis between spin-up 
and spin-down due to vortex metastability. By halting 
the rotation of the normal cloud, we have already found 
that we can reduce the lifetime of the BEC's rotation by a 
factor of 10. Finally, the rotating normal cloud can create 
equilibrated condensates with very large rotation rates, 
which may allow us to approach the regime for which the 
vortices are so close-packed that their separation becomes 
comparable to the healing length j3jj . 
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